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Abstract 

In the field of computer science, the network reliability problem for evaluating the network failure probability 
has been extensively investigated. For a given undirected graph G, the network failure probability is the probability 
that edge failures (i.e., edge erasures) make G unconnected. Edge failures are assumed to occur independently with 
the same probability. The main contributions of the present paper are the upper and lower bounds on the expected 
network failure probability. We herein assume a simple random graph ensemble that is closely related to the Erdos- 
Renyi random graph ensemble. These upper and lower bounds exhibit the typical behavior of the network failure 
probability. The proof is based on the fact that the cut-set space of G is a linear space over F2 spanned by the 
incident matrix of G. The present study shows a close relationship between the ensemble analysis of the expected 
network failure probability and the ensemble analysis of the error detection probability of LDGM codes with column 
weight 2. 



I. Introduction 

A. Background 

Network systems are ubiquitous in communication networks and power line networks, for example. Designing a 
reliable network is critical for achieving a system that is robust to unexpected failures. Theoretical treatments, with 
an abstraction of a real network as an undirected graph, provide insight on robust networks, which will be useful 
for network design. 

The network reliability problem proposed by Moore and Shannon [16] has been extensively investigated primarily 
in the fields of computer science and combinatorics. Assuming that an undirected graph G is given, this undirected 
graph G can be regarded as an abstract model of a real network. For example, in a packet communication scenario, 
a vertex and an edge represent a router and an communication link, respectively. We also assume that each edge 
can be broken. If some edges are broken (the event is referred to as edge failure), then the graph may become 
unconnected. In a packet communication scenario, there exist two routers, between which there is no communication 
route. We refer to such an event as a network failure. Network reliability problems are problems for evaluating the 
probability of a network failure for a given undirected graph. 

Although there are several edge failure models, in the present paper, we adopt the simplest of these. The 
assumption is as follows. Edge failures occur independently, and the probability of an edge failure is uniformly 
given by e(0 < e < 1). We also assume that all of the vertices are perfectly reliable. In the present paper, we focus 
on the all-terminal scenario [7], in which a network is considered to be valid only if there exists a path between 
any two vertices, i.e., the graph representing a network is connected. 

Evaluation of the network failure probability for a given undirected graph is known to be a computationally 
demanding problem. Provan and Ball [17] and Valiant [19] demonstrated that the network reliability problems 
are ^^P-complete, which is a complexity class at least as intractable as MV. Thus, the existence of polynomial- 
time algorithms for solving the network reliability problem appears unlikely. Karger [12] presented a randomized 
polynomial time approximation algorithm for the all terminal network reliability problem. Buzacott [6] and Ball 
and Provan [3], for example, have developed exponential time algorithms for evaluating the exact network failure 
probability with a smaller exponential order. Instead of exact evaluation. Ball and Provan [2] and Colbourn [8] 
presented upper bounds and lower bounds for all-terminal network reliability, which can be evaluated in linear time. 
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B. Main contributions 

In the present paper, we present a probabilistic analysis of all-terminal network reliability problems that is based 
on a random graph ensemble. The random graph ensemble assumed herein is closely related to the Erdos-Renyi 
random graph model [9]. The primary contribution of the present paper is the upper and lower bounds on the 
expected network failure probability. These upper and lower bounds reveal the typical behavior of the network 
failure probability. The upper bound to be presented in Theorem [T] has the following form: 

where k,n are the number of vertices and edges, respectively. The probability Pf{G,e) is the network failure 
probability of G under the condition where the edge failure probability is e. The expectation is taken over a random 
graph ensemble of unconnected graphs with A;- vertices and n-edges. As far as the authors know, this type of bounds 
on the expected network failure probability appear to be novel. 

In order to derive these bounds, an average cut-set weight distribution is derived. The set of cut-set vectors (i.e., 
the incidence vector of a cut-set) of an undirected graph forms a linear- space over binary field F2 spanned by 
an incidence matrix. A combinatorial argument similar to the ensemble analysis of low-density generator matrix 
(LDGM) code with column weight 2 is exploited to obtain the average cut-set weight distribution. The coding 
theory perspective presented in the paper could be useful not only for the network reliability problem but also for 
several graph problems related to the cut-set structure of a graph. 



II. Preliminaries 

A. Random graph ensemble 

We herein consider an ensemble (i.e., a probability space) of undirected graphs. Let k{k > 1) be the number of 
labeled vertices, and let n(l < n < k{k — l)/2) be the number of labeled edges. The vertices are labeled from 1 to 
k, and the edges are labeled from 1 to n. For any adjacent vertices, only a single edge is allowed. In the following, 
[a, b] denotes a set of consecutive integers from a to b. The set Gk,n denotes the set of all undirected weighted 
graphs with A;-vertices and n-edges satisfying the above assumption. For any G G Gfc,n> the sets of vertices and 

edges are denoted by V{G) and E{G), respectively. The cardinality of Gk,n is given by \Gk,n\ = '^•(y)- Here, 
we assign the equal probability 



P{G) = 




for G G „ as the probability measure. Note that the pair {Gk,n,P) defines an ensemble of random graphs and 
is denoted by Gk,n- 



B. Cut-set 

For a given undirected graph G, a partition V{G) = Vi U V2 satisfying = Vi n 1^2 is referred to as a cut of G. 
The set of edges bridging Vi and V2 is referred to as the cut-set corresponding to a cut (Fi, V2). The weight of a 
cut (or a cut-set) is defined as the cardinality of the cut-set and is denoted by uj{E), where is a cut-set. 

For G G Gfc the incidence matrix of G, denoted by M{G) G {0, l}'^'^", is defined as follows: 

, „^ f 1, if ith vertex connects to 7th edge 

^'^^^-^ = I 0; otherwise, ^2) 

where M{G)ij is the (i, j)-element of M(G). 

There are two important properties of M{G). If G is connected, then the rank (over F2) of M{G) h k — 1. The 
row space (over F2) of the incidence matrix M{G) corresponds to the set of all possible cut-sets of G, which is 
called a cut space [5] [10]. Note that the cut space property plays a crucial role in the following analysis. The matrix 
M{G) can be regarded as a generator matrix of LDGM code with column weight 2 [18]. The weight distribution 
of such an LDGM code can be interpreted as the cut-set weight distribution. 
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C. Cut-set weight distribution 

We define the cut-set weight distribution of G by 

B^{G) = ^ l[E is a cut-set of G, oj{E) = v] (3) 

E(ZE{G) 

for non-negative integer v. The function ![•] is an indicator function that takes a value of 1 if the condition is true, 
otherwise this function takes a value of 0. 

The integer-valued function Au^v '■ Gk,n — )■ Z is defined by 

AuAG)= nrnM{G) = c]. (4) 

The set of constant weight binary vectors Z^"''^^ is defined as Z^""'^^ = {x G {0, 1}" : wh{x) = b}, where the 
function wh{-) represents the Hamming weight. The following lemma is another representation of the cut-set weight 
distribution using the incidence matrix of G: 

Lemma 1: Assume G G Gk^n- If G is a connected graph, By{G) is given by 



k 



u=Q 

for u e [0,k],v £ [0,n]. 

Proof: We denote the row space of the incident matrix of G as span(M(G)). The weight distribution By{G) 
can be rewritten as follows: 

B^{G) = l[E is a cut-set of G, u:{E) = v] 

ECE(G) 

Y 1[wh{c) = v]. (6) 

cespan(A'/(G)) 

The second equality is due to the fact that the row space of M{G) is equal to the set of all possible cut-set vectors 
of G. Note that, based on the assumption, the rank of M(G) should he k — 1. The equality 



E ^(^) = ^ E E Hc)I[c = mMiG)] 

c6span(Af(G)) mG{0,l}'' ce{0,l}" 



(V) 



holds for any real-valued function h : {0, 1}" — )■ M because 

|{m G {0,1}'= I mM(G) = c}| = 2 (8) 
holds for any c G span(M(G)). Substituting ^ into (j6]l, we obtain 

MG) = ^ E E nc = mM{G),WH{c) = v] 
me{o,i}'= ce{o,i}" 



k 

^E E E I[c = mM(G)] 

M=OmeZ(''-'") c 
k 

IYAuAG). 



k 



2 

u=0 



(9) 



Note that A„ t,(G) can be seen as an input-output weight distribution of an LDGM code [11] if we regard M{G) 
as a generator matrix with column weight 2. The variables u and v correspond to the input and output weights, 
respectively. Note that Au,v{G) also has a close relationship to the coset weight distribution of the low-density 
parity check (LDPC) code corresponding to a sparse parity check matrix with column weight 2 [20]. 
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III. Ensemble average of cut-set weight distribution 

In this section, we discuss the average (i.e. expectation) of Au^v{G) over the ensemble Gk,n- In the following, 
the expectation operator E is defined as 

E[/(G')]= Yl PiG)f{G), (10) 

where / is any real-valued function defined on Gfc „. 

As a preparation for deriving E[^„ t,(G)], we introduce the following lemma: 

Lemma 2: Assume that m* G Z^'^'") and c* G where u G [0, /c] and v G [0, n]. The following equaUty 

holds: 

E[lKM(G) = c1] = 1 Mfc - n)\ /(^) - - n) Y ^^^^ 



n 



Proof: Due to the symmetry of the ensemble, we can assume without loss of generality that the first u elements 
of m* are one and the remaining elements are zero. In a similar manner, c* is assumed to be a binary vector such 
that the first v elements are one and the remaining elements are zero. 

In the following, we count the number of labeled graphs that satisfy m*M{G) = c* by counting the number of 
binary incidence matrices satisfying the above condition. Let M{G) = (/i /2 • • • fn) where fi is the zth column 
vector of M{G). Since M{G) is an incidence matrix, the column weight of fi is wnifi) = 2 for i G [l,n]. From 
the above assumptions, we obtain 

= ( 1 (12) 

\ 0, i£ [v + l,n]. ^ ^ 



We then count the number of allowable combinations of (/i, /2, • • • , fn) that satisfy ( 12 1. Let 

yl = {/G{0,l}'=|m*/ = l,^x;^^(/) = 2}. (13) 

The cardinality of A is given by |^| = u{k — u) because a non-zero component of / must have an index within 
[1, u], and another non-zero component has an index in the range [u + l,k]. This observation leads to the number 
of possibilities for (/i, /2, . . . , /^), which is given by vl{''^''-''^). The remaining n - v columns, {fv+i, ■ ■ ■ , fn), 
should be taken from the set {/ G {0, 1}*^ | wnif) = 2}\A. Thus, the number of possibilities for such choices is 

summary, the number of allowable combinations of (/i, /2, . . . , /„) is given by 

Y: IKM(G) = c*] =vlin- - ( - - . (14) 



V I \ n — V 
Thus, the left-hand side of ( [TT] ) can be rewritten as follows: 

E[I[m*M(G) = c*] = Y P{G)l[m*M{G) = c*] 

= E lKM(G) = c*] 

v\(n — v)l f u{k — u)^ f (2) — u{k — u] 



V / \ n — V 



(15) 



The final equality is due to ( 14 1. 

The primary result in this section is given as follows: 
Lemma 3: The expectation of Au^v{G) is given by 



where n G [0,k],v e [0, n]. 
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Proof: The expectation of Au^y{G) can be simplified as follows: 

EK,,(G)] = E 



J2 l[mM{G) = c] 
5] J] E[I[mMiG) = c]] 



uj \v 



E [I[m*M(G) = c*]\ , 



(17) 



where the final equality is due to the symmetry of the ensemble. The binary vectors m* G Z^'^'") c* G {0, 1}" can 



be chosen arbitrarily. Substituting ( [TT| ) in the previous Lemma into ( |26| ), we obtain the claim of this lemma. ■ 

IV. Probability of unconnected graphs 

In this section, we discuss the probability such that a randomly chosen graph in Gk^n is unconnected. Such a 
probability has been investigated in detail [9] [4] . These bounds are required for deriving the primary results of the 
present paper, which are described later. 

A. Upper bound on probability for unconnected graphs 

The upper bound on the unconnected probability presented in this section is derived based on the average cut-set 
weight distribution. The unconnected probability Pjj{k,n) is defined as Pjj{k,n) = E [I[G: unconnected]]. An 
upper bound on the unconnected probability is given by the following lemma: 

Lemma 4: The probability for selecting an unconnected graph from the ensemble n is upper bounded by 



k\fC,)-u{k-u) 



u 



n 



1. 



(18) 



Proof: Let pi be the probability such that the incidence matrix of a randomly chosen graph has rank k — i{i G 
[1,A:-1]), i.e., 

Pi= P{G)I[vank{M{G)) = k-i]. (19) 

Using Pi, the unconnected probability Pjj{k,n) can be rewritten as follows: 

Pu{k,n) = E [I [G: unconnected]] 

= Y P{G)I[r3nk{M{G)) < k - I] = ^Pi- (20) 

GeGfc,„ i=2 

The cardinality of the set {m G {0, l}'^ | mM{G) = 0} is denoted by T{G), which can be transformed into 

T{G) = |{m G {0,1}^ I mM(G) = 0}| 

k 

= Y nrnM{G) = 0]. (21) 

«=OmeZ<'''") 

The expectation of T{G) can be expressed simply as follows: 



E[T(G)] 



Y E nmM{G) = 0] 

.M=OmeZ<'=-") 



Y E E[I[mM(G) 

M=Om,eZ('=.") 



E 

u=0 
k 

E 

u=0 



E [I[mM{G) 
(''^''"'^i''^""'') 



(22) 
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The final equality is due to a special case of Lemma |2] The expectation E[T{G)] can be lower bounded by 
2 + 2Pu{k,n) as follows: 

E[T{G)] = 2pi+4p2 + ---2'=-V-i 
> 2pi + 4{p2 + ■ ■ ■ Pk-i) 
= 2{l -Pu{k,n)) + APu{k,n) 

= 2 + 2Pu{k,n). (23) 
This lower bound on E[T{G)] leads to an upper bound on Pu{k,n), as shown below: 

2(y)6^oWV n J 

m 

Note that Erdos and Renyi [9] derived a similar bound based on another combinatorial argument. 

B. Lower bound on probability for unconnected graphs 

The lower bound on the probability for an unconnected graph presented below can be derived based on a simple 
combinatorial argument. 

Lemma 5: The probability of selecting an unconnected graph from the ensemble „ is lower bounded by 

^-('•"'a|(|y(CTO-<'=-K*T'))- 

Proof: The main concept of the lower bound is to bound Pi/{k,n) from below by the probability such that 
a randomly chosen graph G has a zero row vector in its incidence matrix M{G), i.e., 

Pu{k,n) = E [I [G: unconnected]] 

> E [I[M{G) has a zero row vector]] . (26) 

This inequality is obtained for the following reason. If M(G) contains a zero row vector, then G is unconnected 
because G contains an isolated node corresponding to zero vectors of M{G). This implies that 

1[G: unconnected] > 1[M{G) has a zero row vector] (27) 

holds for any G G Gk,n- The right-hand side of ( [26] ) can be further simplified as follows: 

E ^[M{G) has a zero row vector]] = ^ P{G)l[M{G) has a zero row vector 

= — I[M(G) has a zero row vector] 
n!(y) 

> (28) 





n!((P) 

The final equality is due to the following combinatorial argument. Here, we count the number of possible binary 
matrices having a zero vector. Every column of M{G) must contain two-ones, and there are /c — 1-possible positions 
for such two-ones because we assumed that a row are constrained to be zero. This implies that the size of the set 
of possible column vectors becomes ('^2 *) . Therefore, the number of binary matrices having a zero row vector can 
be lower bounded by 

^ l[M{G) has a zero row vector] > n!A; ^(^^ 2 )^ _ _ 1) 2 _ (29) 

The factors k and n\ in the above expression are the number of possible positions of the zero vector and the 
number of possible way to sort n-column vectors, respectively. The negative term —{k — ^ compensates 
the overcounts for multiple zero vectors. Again, this negative compensation term overcounts multiple zero vectors. 



Therefore, RHS of (29) is smaller than LHS of (29). The claim of the lemma follows directly from the above 
discussion. ■ 
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V. Bounds on expected network failure probability 

We assume that the edge failures occur independently with probability e(0 < e < 1). For a given graph G S Gfc,n> 
the edge failures transform G into G' , where G' is referred to as a survivor subgraph of G. The network failure 
probability Pf{G,e) is the probability such that the survivor subgraph G' becomes unconnected. The precise 
definition of the network failure probability is given as follows: 

Definition 1: For G G Gk,n, the network failure probability Pf{G,e) is defined by 

PfiG, e) = I[G: connected] | I[{V{G), E\E'): unconnected]el^'l(l - e)""!-^'! j + I[G: unconnected]. 

\e'ce{g) J 

(30) 

Based on this definition, it is evident that Pf{G,e) = 1 holds if G is unconnected. If G is a connected graph, 
then Pf{G,e) is equal to the probability corresponding to unconnected survivor subgraphs of G. In this section, 
the upper and lower bounds on the expected network failure probability are presented. 



A. Upper bound on expected network failure probability 

The following theorem gives an upper bound on the expected network failure probability. 
Theorem 1 (Upper bound): The expectation of the network failure probability Pf{G,e) over Q^^n can be upper 
bounded by 

2{9)irr,tswv " J\ J 2{fi)).tswv » J 

(31) 



k9) 

Proof: We first derive an upper bound on Pf{G, e) for any G G Gk,n- 



PfiG,e) = I [G: connected] ^ ^\^'): unconnected]el^'l(l - e)""l^'l + unconnected] 

\e'cE{G) J 

< I[G: connected] [ ^ ^ I[E' is a cut-set of G, uj{E') = v^ie + (1 - e))""'' ] + I[G: unconnected] 

\v=l E'CE{G) J 

= I[G: connected] ^^S„(G)e''^ + I[G: unconnected] 

n k \ 

= l[G: connected] - ^ ^ yl„,^(G)e'' +l[G: unconnected] 



2 



^ n k 

- 9 ^u,v{Gy + ^[G- unconnected] . (32) 
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Lemma [T] is used in this derivation. The final inequality is due to the inequality 1[G: connected] < 1 for any 
G G Gfc „. Using the upper bound shown above, the expectation of Pf{G, e) can be upper bounded as follows: 

n k 



E[Pf{G,e)] < E 



2 

D=l lt=0 

^ n k 

9EEeK,-(G)] e- + Pu{k,n) 

,=l u=0 

1 V V (''^ M'' - ""A (^2) - <k -u)\ _^ A (k\ (il) - u{k -u)\_ 
2((P)it.^oWV V )[ n-v +2((P)„eoWV n J"" 

(33) 

The first inequality is due to (32i. Lemmas [3] and |4] are exploited in the derivation. ■ 



+ E[I[G: unconnected]] 



2 

v=l u=0 

n k /,\//, \\ / /k\ /I \\ 1 k 

i - 

< 
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B. Lower bound on an expected network failure probability 

The following theorem includes a lower bound on the expected network failure probability. 
Theorem 2 (Lower bound): The expectation of the network failure probability Pf{G,e) over Qk^n can be lower 
bounded by 

n ^ .... ...X .... ...X , 

n J \ n 



(34) 

Proof: As in the case of the upper bound, we start from the definition of Pf{G, e) and then bound Pf{G, e) 
from below: 

P/(G,e) = I [G: connected] I J] ^\^'): unconnected]el^'l(l - e)""!-^'! | + unconnected] 

\e'ce{g) J 

> I[G: connected] | ^ I[E' is a cut-set of G, uj{E') = v^il - e)''-^ j + I[G: unconnected] 

\i'=l E'CE{G) J 

= I[G: connected] ^^B„(G)e''(l - e)""''^ + I[G: unconnected] 

= I{G: connected] (^^Yl ^«,^(G')e^(l - e)"""] + I[G: unconnected] 

\ v=l u=0 / 

/ n k \ 

= (1 - I[G: unconnected]) U ^u,v{Gy{l - e)"^^ + I[G: unconnected] 

\ v=l u=0 / 
/ n k \ 

- 2 Z Z Au,v{G)e''{l - e)"-^ - -l[G: unconnected] + l[G: unconnected] 

\ ^=1 u=0 / 

(n k \ 
\Y.il ^uAGVil - e)"-" J + h[G: unconnected]. (35) 
v=l u=0 / 

Lemma [T] is used in the above derivation. The inequality based on the binomial theorem 

n k 1 " / \ 1 

v=l u=0 v=0 ^ ^ 

is also exploited in the above transformation. By taking the expectation, we immediately obtain the lower bound 
as follows: 

' ^ n k 

E[Pf{G,e)] > E 



> 



^ n K 

2^Y1 ^u,v{G)e''{l - eY'-^ + -I[G: unconnected] 

^ n k ^ 

- E [A^AG]] e^{l - ey-^ + -Pu{K n) 

(k\ fu{k - u)\ ( (2) - u{k - u)\ _ 



2{<P) 

For bounding Pu{k,n) from below, Lemma |5] has been used 



+ -L((V^M-(A^-i)(V^Ml- (36) 
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VI. Numerical evaluation 



In order to verify the tightness of the bounds on the expected network failure probability proved in the previous 
section, we have performed numerical evaluations for the expected network failure probability. Figure [T] presents 
the upper and lower bounds on the expected network failure probability for the case in which k = 6 and n = 12. 
For the purpose of comparison, the exact values are also plotted in Fig. [T] These exact values have been obtained 
by generating all possible graphs and by evaluating the network failure probability using a recursive computation 
referred to as pivotal decomposition [14]. The horizontal axis represents the edge failure probability e, and the 
vertical axis denotes the value of the expected network failure probability. Note that the upper and lower bounds 
are reasonably tight in this case. Both bounds are very close to the exact value when e is smaller than 10^^. 

Figure [2] shows the values of upper and lower bounds on the expected network failure probability for the case 
in which k = 7 and n = 12. The curve for the exact probability (solid line) approaches a constant value around 
10^^ as e becomes small. The reason for this floor phenomenon is as follows. The unconnected probability is 
Pu{7, 12) = 0.0108 in this case. If e is small, e.g., 10"'^, the dominant failure event is to select an unconnected 
graph from the ensemble. Therefore, the expected network failure probability approaches Pu{7, 12). Note that, in 
the case of Fig. [TJ we can prove Pu{G, 12) = by showing that the upper bound (24i takes a value of zero. Thus, 
we can observe that the curves in Fig. [T] decrease monotonically as e decreases. 
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Fig. I. Expected network failure probability (fc = 6,n = 12): exact, upper and lower bounds 




Fig. 2. Expected network failure probability {k — 7,n = 12): exact, upper and lower bounds 
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VII. Conclusion 

In the present paper, upper and lower bounds on the expected network failure probability are derived. The average 
cut-set weight distribution is key in deriving these bounds. The ensemble analysis used here is similar to the analysis 
of the input-output weight distribution of an LDGM code with column weight 2. The present study reveals a close 
relationship between the ensemble analysis of the network failure probability and the ensemble analysis of the 
error detection probability [21] of LDGM codes with column weight 2. This link between the network reliability 
problem and coding theory may provide a new perspective on the network reliability problem. 

In the present study, we focused solely on the upper and lower bounds for fixed n and k. It would be interesting 
to investigate the asymptotic behavior of the network failure probability when n and k approach infinity while 
maintaining the relationship n = f{k) (f is a real-valued function, e.g., n = fik"^). The error exponent analysis for 
the undetected error probability of LDPC codes shown in [21] is a possible choice for such an asymptotic analysis. 

In the present paper, we discussed a simple graph ensemble, which is closely related to the Erdos-Renyi random 
graph ensemble [9]. In such a graph ensemble, the degree of a node is strongly concentrated around its expectation 
when n and k become large. It seems natural to consider the behavior of the network failure probability for a graph 
ensemble with a non-uniform degree distribution. The analysis for deriving the average cut-set weight distribution 
is approximately equivalent to the analysis of the average weight distribution of an LDGM code ensemble [11] or 
of the average coset weight distribution of an LDPC code ensemble. Some known results, such as those reported in 
[20], may be exploited for the analysis of an ensemble with non-uniform degree distribution. The average cut-set 
weight distribution derived in the present paper may also be useful for several applications other than the evaluation 
of the expected network failure probability. For example, statistical information on the max-fiow of a randomly 
chosen graph could be obtained from the average cut-set weight distribution in combination with the min-cut 
max-flow theorem. 
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